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ABSTRACT 



We investigate the metallicity of the broad line region (BLR) of a sample of 30 quasars in the redshift range 4<z<6.4, by using 
near-IR and optical spectra. We focus on the ratio of the broad lines (SiIV1397+OIV]1402)/CIV1549, which is a good metallicity 
tracer of the BLR. We find that the metallicity of the BLR is very high even in QSOs at z~6. The inferred metallicity of the BLR gas 
is so high (several times solar) that metal ejection or mixing with lower metallicity gas in the host galaxy is required to match the 
metallicities observed in local massive galaxies. On average, the observed metallicity changes neither among quasars in the observed 
redshift range 4<z<6.4, nor when compared with quasars at lower redshifts. We show that the apparent lack of metallicity evolution 
is a likely consequence of both the black hole-galaxy co-evolution and of selection effects. The data also suggest a lack of evolution 
in the carbon abundance, even among z>6 quasars. The latter result is puzzling, since the minimum enrichment timescale of carbon 
is about 1 Gyr, i.e. longer than the age of the universe at z~6. 
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1. Introduction 

One of the important results from studies of high-redshift 
quasars is the apparent lack of evolution in their emission proper- 
ties. In particular, the lack of evolution in the emission-line ratios 
is interpreted as a lack of metallicity evolution among QSOs at 
different redshifts. Previous tentative indications of metallicity 
evolution with redshift were later ascribed to a metalli city de- 
pendence on the quasar luminosity (Naga o et al.Ll2006ah . which 
may mimic a dependen ce on redshift i n flux-limited samples. 
IWarner et all (2004) and lShemmer et all ([2004) suggest that the 
metallicity-luminosity relation may actually result from a BH- 
metallicity dependence or from a relation between accretion rate 
(L/Lndd) and metallicity. 

At z>5 most of the emission lines used to constrain the BLR 
metallicity are redshifted into the near-IR. As a consequence 
the investigation of the BLR metallicity at z>5 has been lim- 
ited to a small number of objects. Various authors have investi- 
gated the iron abundance in z~6 QSOs by using near-IR spec- 
tra to measure the intensity of the UV "Fell-bump" (Arest ~ 



and NV1240 emission lines, suggesting super-solar metallici- 
ties. However, the calibration of the ratio NV/CIV as a tracer 
of gas metallicity is mostly based on the assumption that nitro- 
gen is a secondary element and that, therefore, its abundance 
scales quadratically with metallicity, i.e. N/H oc (O/H)^. The lat- 
ter assumption may, however, be an oversimplification and may 
give unrealistically high metallicities, since the nitrogen abun- 
dance actually depends on the detailed star formation history 
jBresoUn et aUl2004tljiang et al.Ll2008l) . 



2000 - 3050A) relative to flie Mgm2789 line dFreudUng et al 



2003'; ' Maiolino et all l2003t llwamuro et all |2004 iKurk etal 



2007), as a proxy of the Fe/a-elements abundance ratio. No evo- 
lution was found in the ratio Fell/Mgll; however, since the Fell 
bump is a strong BLR coolant, its strength depends on several 
other parameters and may not depend primarily on the iron abun- 

darice. 

iPentericci et al.l ( |2002|) investigated the metallicity in the 
BLR of two QSOs at z~6 by measuring the ratio of the CIV 1549 



* Based on data obtained at the VLT, through the ESO programs 
074.B-0185 and 075.B-0487, and at the Telescopio Nazionale Gahleo. 



iNagao et al] (l2006ab show that the ratio 
(SiIV/ll397H-OIV]^1402)/CIV/ll549 (hereafter 
(SiIVH-OIV)/CIV) is a tracer of the BLR metallicity, the 
reason mainly being that the relative importan ce of CIV as 
a coo lant decreases with the BLR metallicity iFerland et al.l 
( Il996l) . This method has the advantage of not relying on the 
assumption N/H oc (O/Uf. Moreover, (SilV-i-OIV) and CIV are 
well isolated features, which do not require difficult deblending 
techniques a nd which can be measured even in low-resolution 
spectra. Jian g et aT] (l2007 h used near-IR spectra to measure 
(SiIV-i-OIV)/CIV in five QSOs at z~6, and find no evidence of 
evolution when compared with lower redshift QSOs. 



In this paper we present measurements of the 
(SiIV-HOIV)/CIV ratio in 30 QSOs at 4<z<6.4 obtained 
by means of near-IR spectra, with the goal of better investigat- 
ing the metallicity evolution of QSOs in the early Universe. In 
the following we adopt the following cosmological parameters: 
//o=70 km S-' Mpc-', nM=0.3, and 
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2. Observations 

We observed a sample of 30 high-redshift quasars (4.0 < 
z < 6.4) from the SDSS by means of near-IR and opti- 
cal spectra covering at least the UV rest-frame emission lines 
SiIV^1397+OIV]/tl402 and CIV^1549, but in most cases the 
spectra extend to /Irest ~ 3000 - 4000A. The original goal of 
most of the observations was to constrain the dust extinction in 
high-z QSOs. A more detailed description of the data and the 
results on the dust extinction will be given in Gallerani et al. (in 
prep.). Here we only focus on a byproduct, namely the evolution 
of the BLR metallicity based on the (SiIV+OIV)/CIV ratio. 

Observations were obtained both with the Italian Telescopio 
Nazionale Galileo (TNG) in Spain and with the Very Large 
Telescope (VLT)-ESO in Chile. Observations were performed 
in several observing runs from 2003 to 2005. The observa- 
tions at the TNG were obtained with the Near Infrared Camera 
Spectrograph (NICS) mostly with the Amici prism to obtain 
spectra in the range 0.9-2.3 jum at R~75. This low-resolution 
mode is excellent for investigating the QSO continuum shape, 
but also for detecting broad emission lines. Some QSOs were 
observed again with the IJ grism to obtain 0.9- L45 yum spectra 
at R~500. Typical integration times ranged from ~20 minutes to 
~3 hours. The obser ving strategy and data reduction are similar 
to those discussed in lMaiolino et al.l (|2004|) . 

The spectroscopic observations at ESO-VLT were done with 
the FORS2, along with the grism GRIS 1501, to observe the range 
6000- 11 000 A at R~300. These observations are mostly used to 
cover the short-wavelength part of some of the quasar spectra 
not properly sampled by the near-IR observations, but we also 
specifically observed a few quasars with no near-IR data with 
the specific aim of measuring the (SiIV+OIV)/CIV ratio. The 
total exposure times range from 30 to 60 minutes. For some of 
the z < 5 quasars observed with NICS, for which no FORS2 ob- 
servations were availab le, we combine our ne ar-IR spectra with 
optical data taken from I Anderson et al] (1200 lb . 
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Fig. 1. {Upper) (SiIV-i-OIV)/CIV flux ratio as a function of red- 
shift for quasars observed in this work. Squares indicate non- 
BAL quasars, while diamonds indicate BAL quasars. {Lower) 
(SiIV+OIV)/CIV ratio inferred from quasar stacked spectra, 
both for the high-redshift quasars presented here (filled symbols) 
and for the lower redshift quasars studied in lNagao et al.l (l2006ah 
(empty symbols). 



3. Results 

While individual spectra will be shown in a forthcoming 
paper (Gallerani et al., in prep.), in this letter we report 
the measurements of the (SiIV-FOIV)/CIV ratio. To mea- 
sure emission line fluxes, these low resolution quasar spectra 
were fitted with a model consisting of a power-law contin- 
uum and Gaussian profiles to model the emission lines. The 
power-law continuum was fitted using small spectral regions 
free of emission features (1345-1355A,1445-1455A,1695- 
1704A,1973-1983A). The SilV+OIV] blend is unresolved and 
is fitted with a single Gaussian. Eight of our sources are broad 
absorption line (BAL) quasars. In these cases the measurement 
of the emission line flux is more complex and uncertain be- 
cause of the absorption troughs, which are avoided when fit- 
ting both the continuum and the emission lines. Table [T] lists 
the (SiIV-i-OIV)/CIV ratio for all quasars in our sample, along 
with their redshift and continuum luminosity {AL^ at 1450A). 
Figure [1] (top panel) shows the measured (SiIV+OIV)/CIV ra- 
tios as a function of re dshift. According to the calibration given 
in Nagao et alj (l2006al) . the observed ratios correspond to metal- 
licities that are several times solar. No correlation is observed 
between the emission line ratio (hence metallicity) and red- 
shift. Figure|2]shows the average spectra in the redshift intervals 
4.0<z<5.5 and 5.5<z<6.4 (obtained by first normalizing indi- 
vidual spectra to the continuum at /11450A, and avoiding BAL 
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Fig. 2. Stacked spectra of quasars in different redshift bins (all 
normalized to A - 1450A). The red solid line and the blue 
dashed line are the stacked spectra from our sample in the red- 
shift intervals 5.5<z<6.4 and 4.0<z<5.5, respectively. The black 
dotted line is the SDSS quasars stacked spectrum in the redshift 
interval 2.5<z<3.5. 

quasars), which are remarkably similafl again highlighting any 
lack of evolution. 

' Small differences in Lye and NV are due to different absorption of 
the IGM and to different spectral resolution of the spectra. 
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Table 1. (SiIV+OIV)/CIV measurements. 



Name 



(SiIV+OIV)/CIV log^L;', 



SDSS J000239.4+255035 


5.80 


0.61+0.12 


46.93 


SDSS J000552.3-000656 


5.85 


0.34+0.06 


46.20 


SDSS J001714.66-100055.4'' 


5.01 


0.58+0.07 


46.56 


SDSS J012004.82+ 141 108.2'' 


4.73 


0.32+0.03 


46.09 


SDSS J015642.1 1+141944.3'' 


4.32 


0.16+0.03 


46.62 


SDSS J023 137.6-072855 


5.41 


0.29+0.07 


46.56 


SDSS J023923. 47-08 1005.1'' 


4.02 


0.23+0.04 


46.58 


SDSS J033829.3+002156 


5.00 


0.97+0.29 


46.46 


SDSS J075618.1+410408 


5.07 


0.80+0.15 


46.53 


SDSS J083643.8+005453 


5.80 


0.63+0.16 


46.97 


SDSS J085210.89+535948.9'' 


4.22 


0.42+0.04 


46.53 


SDSS J095707.67+061059.5 


5.16 


0.70+0.14 


46.65 


SDSSpJ1021 19.16-030937.2 


4.70 


0.35+0.12 


46.58 


SDSS J103027. 1+052455 


6.28 


0.59+0.20 


46.68 


SDSS J104433.04-012502.2'' 


5.78 


0.40+0.13 


46.88 


SDSS J104845.05+463718.3'' 


6.20 


0.42+0.20 


46.81 


SDSS J114816.6+525150 


6.40 


0.41+0.08 


46.95 


SDSS J120441.7-002150 


5.05 


0.57+0.12 


46.63 


SDSSpJ120823.8+001028 


5.27 


0.63+0.21 


46.07 


SDSS J130608.2+035626 


5.99 


0.38+0.19 


47.32 


SDSS J141111. 3+121737 


5.93 


0.57+0.13 


46.58 


SDSS J160254.2+422823 


6.07 


0.53+0.17 


46.90 


SDSS J160320.89+072104.5 


4.39 


0.39+0.08 


46.85 


SDSS J160501.21-01 1220.6'' 


4.92 


0.37+0.15 


46.46 


SDSS J161425. 13+464028.9 


5.31 


0.34+0.03 


46.62 


SDSS J162331. 8+31 1201 


6.22 


0.54+0.12 


46.54 


SDSS J162626.50+275132.4 


5.20 


0.30+0.09 


46.94 


SDSS J163033.9+401210 


6.06 


0.23+0.11 


46.38 


SDSS J220008.7+001744 


4.77 


0.55+0.07 


46.58 


SDSS J221644.0+001348 


4.99 


0.29+0.06 


46.18 



" Rest-frame continuum luminosity at 1450 A, in units of erg s 
'' BAL quasars. 



Our results can be compared with the line ratios and metalhc- 
ities inf erred at even lower redshifts, as obtained by Nagao et al.l 
(l2006al) . It is important to use quasars with similar luminosi- 
ties, since quasar metallicities show a strong correlation with 
luminosity. Given that flux limited surveys probe different lu- 
minosities at different redshifts, this may introduce an appar- 
ent redshift evolutio n. Therefore, from the 2<z<4 templates in 
iNagao et al.l (l2006al) . we select those in the highest luminosity 
bin (-28.5 > Mb > -29.5), which match the quasar lumi- 
nosities in our sample (which are also restricted to a relatively 
narrow range). In Fig. |2] the average spectra of z>4 quasars 
are compared with the average spectrum of luminous QSOs at 
2.5<z<3.5, which are all very similar The lack of any evolu- 
tion with redshift can be also appreciated in the bottom panel of 
FiglT] which shows the ratio (SiIV+OIV)/CIV measured in aver- 
age quasar spectra at different redshifts (matched in luminosity). 
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Fig. 3. Evolution of the gas metallicity (top), bolometric quasar 
luminosity (middle) and total gas mass (bottom) in the model by 
Granato et al.. (2004) . The horizontal dashed line shows the ap- 
proximate minimum bolometric luminosity detectable at high-z 
in the SDSS quasar survey. The shaded area shows the epoch 
when a high-z quasar is luminous enough to be detected in the 
SDSS survey (if unobscured). The hatched area is the epoch 
when the quasar has ejected more than half of the gas mass, 
and it is likely unobscured and therefore detectable in optical 
surveys. By the time the quasar has reached the luminous, "un- 
obscured" epoch the gas is already highly enriched. 



4. Discussion 

4.1. The quasar metallicity "evolution" 

The apparent lack of evolution observed in Figs. [TJ[2]should not 
be interpreted as a lack of evolution of the BLR metallicity in 
individual AGNs. Indeed, Fig. [1] shows the average metallicity 
of the BLR in quasars that are accreting at the given redshift, but 
does not trace the evolutionary path of individual quasars. The 
apparent lack of evolution in the BLR metallicity observed in 
Fig. [1] likely results from a combination of the BH-galaxy co- 
evolution and selection effects. Indeed, to cross the detection 
threshold of the SDSS magnitude-limited survey, high-redshift 



quasars must have high luminosities, hence (even if accreting at 
the Eddington limit) high black hole masses. Most models pre- 
dict that high black hole masses must have been accompanied b 
th e formation o f a massi ve host galaxy (e^. Granato et al., 200' 
iDi Matteo et all 120051: iHopkins et al.[ 120081: iLi et all l200' 
which would result into the local Mbh - Mspheroid relationship. 
Therefore, by the time a quasar at any redshift is detectable in 
a magnitude-limited survey, its host galaxy must have evolved 
significantly and enriched its ISM significantly. The quasar feed- 
back is another evolutionary effect that may yield to observa- 
tional biases resulting in an apparent lack of metallicity evolu- 
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tion. Indeed, according to many models, during the early phases, 
when the host galaxy is still metal poor, the accreting black hole 
is embedded within the dusty ISM, and therefore difficult to de- 
tect in optical surveys. Only during the late evolutionary phases, 
when the galaxy is already metal rich, the quasar develops winds 
powerful enough to expel large quantities of gas and dust, so that 
the quasar becomes visible to optical observations. 

To show these combined effects in a more quantitative way, 
w e exploited the quas ar-galaxy evolutionary models reported 
in iGranato et alj (l2004i) . In Fig. [3] we show the evolution of 
the gas metallicity in a massive galaxy forming in a dark halo 
of Mhaio = lO'-''^ Mq, along with the evolution of the quasar 
bolometric luminosity and the total gas mass. While the host 
galaxy forms stars, the gas metallicity increases. At the same 
time the black hole accretion increases, yielding an increasingly 
high quasar luminosity. At z>5 the SDSS magnitude limit trans- 
lates into a minimum quasar bolometric luminosity of about 
10"*^ erg s"' for detection. Fig. [3] shows that by the time the 
quasar reaches such high luminosity (shaded region), the metal- 
licity in the host galaxy is already higher than ~ 3 Z©. At this 
time the model expects the quasar to be probably still embed- 
ded in gas and dust, although the quasar has already developed a 
wind expelling gas. It is not easy to identify the stage when the 
quasar becomes unobscured and detectable in optical surveys, 
since it depends on the detailed distribution of dust. We assume 
that the quasar becomes optically visible when more than half 
of the gas mass has been expelled. This is a somewhat arbitrary 
assumption, but we not e that i t is roughly consistent with what 
obtained by iLapi et al] (|2006|) . who identify the epoch of un- 
obscured accretion as starting a few 10^ yrs before the epoch 
of maximum luminosity. In Fig. |3]we note that by the time the 
quasar is unobscured (hatched region), the gas in the host galaxy 
has already reached a metallicity of about ~ 4 Z©. Summarizing, 
the co-evolution of black holes and galaxies, combined with ob- 
servational selection effects (mostly in optical surveys), naturally 
explains the finding that unobscured quasars of a given luminos- 
ity appear to have on average the same metallicity at any redshift. 

4.2. The extreme metallicities in the BLR 

According to iNagao et al.l (l2006 al). the ratio 
(SiIV+OIV)/CIV~0.4 observed in the stacked spectrum of 
the most distant quasars corresponds to a gas metallicity 
of ~ 7 Zq . Such huge metallicities were also inferred by 
iNagao et al.l (l2006al) based on a much wider set of lines of lower 
redshift quasars. These high gas metallicities are not unrealistic. 
Indeed, the BLR is a small nuclear region (less than a few pc in 
the most luminous QSOs) with masses of a few times ~ 10"* M©, 
which can be enriched in situ to super-solar metallicities within 
less than 10**yrs by having a supernova explosion less often than 
every 10'*yrs. More troublesome is that such high metallicities 
are not found in the stellar population of local massive galaxies, 
not even in their central regions. Therefore, the high-metallicity 
gas observed in the BLR of high-z quasars must either be 
expelled (as in quasar feedback models) or diluted by lower 
metallicity gas in the host galaxy before forming stars. 

A related issue is how much the BLR metallicity is repre- 
sentative of the gas metallicity in the host galaxy. Investigating 
this issue would require accessing metallicity diagnostics that 
directly trace the quasar host galaxies. Host galaxies have been 
probed either through the narrow emission lines (only observ- 
able in obscured QSOs, e.g. Nagao et al., 200 6b) or through th e 
associated narrow absorption lines (e.g. lD'Odorico et al.Ll2004l) . 
but only at intermediate redshifts (z<3). It was found that the 



metallicity of the host galaxy is typically lower than observed in 
the BLR, but still super-solar, and again without any evidence 
of evolution with redshift. No direct measurements are currently 
available for the metallicity in the host galaxy of z~6 quasars . 
Howe ver, the detection of huge quantities of dust jBeelen et all 
120061) suggests that even the host of these quasars are heavily 
enriched. 

4.3. The carbon problem 

Relative elemental abundances are much more difficult to infer 
However, the observational result that the ratio (SiIV+OIV)/CIV 
does not evolve suggests that the abundance of carbon relative to 
silicon and oxygen also does not evolve significantly. Large car- 
bon abundances in the hosts of the most distant quasars is also 
suggested by the detection of strong [CIII lSSAt m and CO lines 
in the most distant quasars dMaiolino et alil2005tlBertoldi et all 
|2003>) . Carbon is mostly produced by AGB stars and planetary 
nebulae, most of which evolve on long timescales. Although the 
first AGB stars appear as soon as ~50 Myr after the onset of 
star formation, the bulk of carbon production occurs after about 
1 Gyr, yielding a delayed enrichment with respect to a-elements 
(e.g. oxygen) that are promptly produced by SNifl At z>6 the 
age of the universe is less than 1 Gyr, hence stellar evolution may 
fall short of the time needed to produce the observed large car- 
bon abundance. We note that this issue is independent of the size 
and mass of the BLR, making it just a pure timescale problem. 
Tackling this issue requires a more accurate determination of the 
carbon abundance, which may come from future high spectral 
resolution optical/near-IR observations o r from fut ure submm 
observations of far-IR fine structure lines (iMaiolinoli200 8). 
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